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Abstract. We have investigated the spatial distribution, and the properties and chemical composition of the 
dust orbiting HD 100546. This system is remarkably different from other isolated Herbig Ae/Be stars in both 
the strength of the mid-IR excess and the composition of the circumstellar dust. To explain spectral features 
and the amount of mid-IR dust emission the presence of a component of small (< 10 (im) grains radiating at 
~200 K is required, which is not seen in other well investigated Herbig Ae/Be systems. This additional component 
■^j- ' is inconsistent with a uniform flaring disk model. The fraction of intercepted stellar light that is absorbed and 

l/^ , re-emitted in the mid-IR is so large (~ 70%) that it requires the disk to be more 'puffed up' at about 10 AU, 

where the grains have T ~ 200 K. This may occur if a proto- Jupiter clears out a gap at this distance allowing 
direct stellar light to produce an extended rim at the far side of the gap. The other remarkable difference with 
other isolated Herbig Ae/Be systems is the presence of a much larger mass fraction of the crystalline silicate 
forsterite in the circumstellar dust. We find that the mass fraction of crystalline silicates in HD 100546 increases 
with decreasing temperature, i.e. with larger radial distances from the central star. This distribution of crystalline 
dust is inconsistent with radial mixing models where the crystalline silicates are formed by thermal annealing 
l' above the glass temperature in the very inner parts of the disk, and are subsequently transported outwards and 

1~ \ mixed with amorphous material. We speculate that the formation and spatial distribution of the crystalline dust 

may be linked to the formation of a proto- Jupiter in the disk around HD 100546. Such a proto- Jupiter could 
gravitationally stir the disk leading to a collisional cascade of asteroidal sized objects producing small crystalline 
grains, or it could cause shocks by tidal interacti o n wit h the disk which might produce crystalline dust grains 
through flash heating. As shown lw|Ma lfait et al.| <|l998T) . the infrared spectrum of HD 100546 is very similar to 
that of C/1995 Ol Hale-Bopp llCrovisier et all Il997t) . Using an identical methodology, we have therefore also 
studied this solar system comet. Both objects have an almost identical grain composition, but with the important 
difference that the individual dust species in Hale-Bopp are in thermal contact with each other, while this is not 
the case in HD 100546. This suggests that if similar processes leading to the dust composition as seen in HD 100546 
also occurred in our own solar system, that Hale-Bopp formed after the formation of one or more proto-gas giants. 

Key words. Circumstellar matter - Stars: formation - Stars: pre-main-sequence 



1. Introduction 

Herbig Ae/Be stars (hereafter referr ed to as HA EBE 
stars) were first described as a group bv lHerbid (|l960) , in a 
study which was aimed at finding intermediate mass young- 
stars. Many studies have since confirmed the young pre- 
main-sequence (PMS) nature of HAEBE stars. Hipparcos 
parallaxes (e.g. van den Ancker et al. 1998) allowed accu- 
rate positions in the HR diagram to be obtained, and from 
comparison to PMS evolutionary tracks masses of the or- 
der of 2 to 8 M Q and stellar ages measured relative to 
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the birthline of typically 10 6 to 10 7 yrs were found. Direct 
imaging at millimeter wavelengths of several Herbig Ae 
stars revealed the presence of rota ting flattened structures 
l|Mannings fc Sargentll997L EoOOli. believed to be the rem- 
nant of the accretion disk and the site of on-going planet 
formation. It is not clear whether disks are also common 
around more massive Herbig Be stars. 

At infra-red ( IR) wavelengths, t he Infrared Space 
Observatory (ISO iKessler et alJ Il996j) has obtained full 
2-200 rim spectra of the brightest HAEBE stars. These 
spectra show a wealth of detail concerning the thermal 
emission of dust in the circumstellar environment. In a 
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series of paper s, we have studied the dust emission from 
HAEBE stars JWaelkens et al.lll996t IMalfait et alJll998l 
1999[ Ivan den Ancker et alJ 1200(1 l200lt iBouwman et all 
2000ll200lHMeeus et al.l200l[) . focusing on the mineralog- 



ical composition of the dust. These studies indicate that a 
substantial modification of the dust composition from that 
in the interstellar medium (ISM) occurs on time scales 
that are still poorly constrained, but are generally less 
than 10 7 years. We mention grain growth, the formation 
of crystalline silicates and of crystalline H2O ice, and of 
hydro silicates. These modifications are important clues 
to the processes that eventually lead to planet formation, 
and which can be compared to the records of planet forma- 
tion as found in solar system objects (meteorites, comets, 
interplanetary dust particles). 

In this study we re-exa mine the ISO s pectr um of 
HD 100546, first presented bv lWaelkens et alJ (|l996|) . This 
object has an exceptionally high fraction of crystalline 
silicates, and its ISO spectrum shows a remarkable re- 
semblance to that of t he solar system com et Hale-Bopp 
( Crovisier et"ai]ll997() . IMalfait et all l)l99S|) already sug- 
gested the exciting possibility that we are witnessing the 
birth of an Oort cloud with a multitude of cometary 
bodies being scattered by a (hypothesized) giant planet. 
IBouwman et alJ l(200li ) analyzed the 10 /im silicate band 
of 14 HAEBE stars including HD 100546, and found that 
its mineralogical composition deviates substantially from 
that of other HAEBE systems: in addition to a large frac- 
tion of forsterite, the average grain size of the particles 
causing the 10 /mi emission is larger. Also, a lack of sil- 
ica compared to forsterite was found. The latter can be 
interpreted as a difference in the crystallization process 
leading to the formation of forsterite. We decided that 
a closer look at the distribution of the different miner- 
alogical components in terms of mass and temperature is 
needed to better understand the nature of HD 100546. We 
also analyze the ISO spectrum of Hale-Bopp, using the 
same analysis method as for HD 100546. Previously, we 
have carried out a similar detailed analysis for AB Aurigae 
and HD 163296 l(Bouwman et alJl2000h and we will com- 
pare these results with those obtained here for HD 100546. 
The obser vations of HD 10054 6 use d in this paper ar e 
taken from IMalfait et all l)l99S|) and lMeeus et all l|200lj) . 
The observations of Hale-Bopp were secured when the 
comet had a he l iocent ric distance of 2.8 AU and are from 
ICrovisier et all {l997h. 

This paper is organized as follows: in Sect. [21 we dis- 
cuss the difference in dust composition and spectral en- 
ergy distribution of HD 100546 compared to other HAEBE 
systems. In Sect.|3we explain the method to fit the spec- 
trum. Sect. ^describes our results for HD 100546 and for 
Hale-Bopp. In Sect. we compare our results to those ob- 
tained for AB Aur and HD 163296, and discuss the impli- 
cations for the formation process of the crystalline silicates 
in HD 100546. Sect.HOsummarizes the results of our study. 
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Fig. 1. Comparison between the energy distribution of 
AB Aur (light line) and HD 100546 (dark line). The spec- 
tra have been normalized to the respective stellar lumi- 
nosities. Plotted are the combined ISO-SWS and LWS 
spectra. Note the large difference in the near- and mid- 
IR luminosity of both systems. 



2. The deviating dust composition and spectral 
energy distribution of HD 100546 

The spectral appearance and dust properties of HD 100546 
differ from that of other HAEBE stars. We will start out 
with a summary of these differences: 

First, the ISO-spectra of HD 100546 show strong and 
pronounced emission, identified with the crystalline sili- 
cate forsterite, which in abundance m ust be far in excess 
of that seen in oth er HAEBE systems IjMalfait et alJ ll998: 
iMeeus et alJl200ll). 

Second, IBouwman et alJ l|200lf) identified an emission 
component at 8.6 /xm in the spectra of a large sample of 
HAEBE stars which they attribute to silica (Si02), and 
which seems to be correlated with forsterite. This corre- 
lation is interpreted as evidence that the crystalline sili- 
cates are formed by thermal annealing. However, the dust 
properties in HD 100546 - and in this respect also that of 
comet Hale-Bopp - does not comply with this correlation 
between silica and forsterite. This suggests that a different 
formation process is responsible for the large mass fraction 
of forsterite seen in both objects. 

Thir d, from a detailed an alysis of the 10 fim spectral 
region <|Bouwman et alJl200fl . it was shown that the typ- 
ical grain size of the silicate particles emitting at these 
wavelengths is larger compared to the sizes found in other 
HAEBE systems. This suggests the disk in HD 100546 is 
more evolved. 
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AB Aur 
HD 100546 



Ld/L* 
0.46 
0.51 



Lnir/L(J 

0.51 
0.17 



Lmir/L<i 
0.34 
0.70 



LFm/Ld 
0.15 
0.13 



Table 1. Luminosity of the circumstellar dust in 
AB Aur and H 100546 as a fraction of the stellar luminos- 
ity (first column). The remaining columns list fractional 
dust luminosities in the near-IR (A < 8/im), the mid-IR 
(8 < A < 50/xm) and far-IR (A > 50^m), respectively. 



Fourth, the spectral energy distribution of HD 100546 
shows remarkable differences when compared to other 
HAEBE stars. To illustrate the difference, we have plot- 
ted in Fig. ^ the SEDs of AB Aurigae (light line) and 
HD 100546 (dark line) normalized to the stellar luminos- 
ity. Both dust disks intercept about 50% of the stellar 
light (see Tablc^). It is easily seen that ABAur re-emits 
a larger fraction of this intercepted light in the near-IR, 
while HD 100546 shows a larger emission at mid-IR wave- 
lengths. The far-IR and sub-millimetre luminosities are 
comparable. The difference in the way radiation is redis- 
tributed by the dust essentially implies that the stellar 
flux is absorbed at different locations in the circumstellar 
disk. Dust grains dominating the near-IR emission have 
temperatures of at least ~ 1000 K. For a grain to reach 
such a temperature it has to be within ~ 1 AU of the 
central star. In the case of ABAur it is within this region 
that about half of the intercepted light is reprocessed (see 
Tabled . HD 100546 emits a relatively modest amount of 
radiation in the near-IR, however, 70% of the total dust 
luminosity comes out in the mid-IR. Grains emitting at 
these wavelengths have typical temperatures of ~ 200 K. 
Assuming these grains receive direct stellar light this im- 
plies that the re-emitted mid-IR flux originates from about 
10 AU from the central star. 

In estimating the above emitting regions, we have as- 
sumed that the dust medium is optically thin. Let us, 
for the moment, hold on to this assumption. For r <C 1 
the emission scales with dust mass, implying that the 
disk of HD 100546 is less massive in the inner parts 
compared to ABAur. Giv en that the system is older 
Jvan den Ancker et alll99fil> . it suggests that a larger frac- 
tion of the circumstellar matter has been removed from the 
disk (evaporated or accreted) and/or is incorporated into 
larger bodies, substantially reducing the near-IR emission. 
In an optically thick view, the emissio n scales not with 
the du st mass but with emitting surface. fPullemond et alJ 
(2001) show that the near-IR excess may be due to a 
puffed-up inner rim surface. The height of this inner rim is 
a function of the stellar luminosity and of the surface den- 
sity in the disk near the evaporation temperature of the 
dust. The low near-IR excess of HD100546 compared to 
AB Aur is then the combined effect of the smaller stellar 
luminosity, and a severe depletion of the inner disk. 

Co ncerning the mid- I R flux es, similar considerations 
apply. IDullemond et all ((2001) show that the SED of 



AB Aur can be reproduced with a passive reprocessing 
flaring disk in hydrostatic equilibrium. As the mid-IR 
spectrum of t his star is similar to that of many other 
HAEBE stars l(Meeus et al.ll200l|) . most likely this model 
has generic value. However, the model can not explain the 
SED of HD 100546. For a uniform flaring disk model to 
reproduce the strong mid-IR emission, the opening angle 
of the disk has to increase faster with radius in order to in- 
tercept more stellar light relatively c lose to the star. Given 
the comparable sizes of the systems ((Augereau et all200lt 
iMannings fc Sargentlll997j) . such a solution would lead to 
a much larger total covering fraction. However, these are 
about 50% in both systems. 

A possible solution for this problem is to induce a 
sizable gap in the circumstellar disk. Such a gap would 
produce a large rim at the far side of the gap (when 
viewed from the direction of the central star), of which 
the largest part would be irradiated directly by the star 
(see Fig. EH- As a result, a larger absorbing surface is 
created and consequently the larger emitting surface that 
is required. IDullemond et all l|200lh show that in such a 
geometry indeed the disk scale height is increased sub- 
stantially, increasing the emission from the region where 
the dust is ~ 200 K. In Sect. I5.ll we will speculate that a 
massive planet, which has formed in the disk at ~ 10 AU, 
may have created such a gap. 

We conclude that HD 100546 differs considerably from 
other HAEBE systems in both dust composition and disk 
geometry. This implies that processes different from those 
observed in other HAEBE systems currently dominate the 
HD 100546 system. In this paper we will try to trace these 
processes. 

3. Modelling 

The diagnostic method adopted is identical for the Herbig 
Be star HD 100546 and comet Hale-Bopp. We use the ra- 
diative transfer program modust to model the circumstel- 
lar dust using the mode in which the material is assumed 
to be optically thin. This is correct for comet Hale-Bopp, 
but for the circumstellar disk of HD 100546 this may not 
be valid. However, in discussing dust properties we will 
focus on the mass over temperature distribution of the 
material responsible for the infrared emission. This diag- 
nostic does not depend on the optical depth properties of 
the medium and therefore the simple approximation used 
here is still meaningful. For both objects the material is 
distributed in a spherical shell, which for Hale-Bopp is 
positioned at the distance from the sun corresponding to 
that of the location of the coma and tail at the time of 
observation. For recent applications and for descriptions 
of tec hniq ues used in MODUST see e.g. iBouwman et alJ 
l|2000l) and lKemner et al.l l|200lh . 

Th e grains are irr adiated by the central star, for which 
we use iKurucd (fl993) energy distributions. The Hipparcos 
distance to the B9Vne star HD 1 00546 is 103 ± 7 pc, 
yield ing a luminosity L — 32 Lq l|van den Ancker et alJ 
1998). Comparison with evolutionary tracks, using T c g = 
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10 500 K, places the star on the main sequence having an 
age ol > 10 Myr and mass of 2.4 Mq . In modelling Hale- 
Bopp, we use Kurucz's solar model with T e g = 5 777 K. 

The absorption of stellar radiation at mostly ultravio- 
let and optical wavelengths and re-emission in the infrared 
is consistently taken into account assuming the particles 
are in radiative equilibrium, yielding the spectral energy 
distribution. A key difference of this method compared to 
other often employed approaches in modelling circumstel- 
lar and cometary spectra is that both the spectral char- 
acteristics responsible for the distinct features as well as 
those responsible for the featureless continuum are mod- 
elled simultaneously. Other methods usually introduce an 
artificial separation between continuum and features, fit- 
ting the spectrum combining empiric al temperatures with 
laboratory extincti on efficiencies fe.g. lWooden et alJl 999: 
iMason et al.l l200lh or with template extinctions based 
on measurements of inter planetary dust particles (IDP's) 
(e.g. IWooden et ai]l2000lh complicating a reliable inter- 
pretation. 

In the near- and mid-IR, typically carbon and large 
10 yum) amorphous silicate grains are responsible for 
the overall continuum, while water ice is a dominant con- 
tributor in de far-IR and sub-millimetre range. At 10 /im 
the molecular vibrational modes of Si-0 bonds in small 
amorphous silicate grains 10 fim) produce a distinct 
feature, as do silicate crystals in the mid-IR part of the 
spectrum. These latter particles, however, also produce 
significant flux outside of their resonances, which one can 
not a priori distinguish from continuum only contributors. 
Continuum subtraction methods therefore are susceptible 
to (systematic) errors. This may significantly affect tem- 
peratures, abundances, and sizes derived for especially the 
crystalline material. Our method does not suffer from this 
problem. In fairness, we should mention that in other re- 
spects studies using ground based data have advantages 
over the interpretation of ISO data. For instance, the mid- 
IR observations analysed by Wooden al. (2000) - which 
have been obtained close in time to the ISO data - have 
much smaller beams, allowing one to study gradients in 
the particle properties as a function of distance to the 
cometary nucleus. 

3.1. Size and shape properties of grains 

Composition, size and shape properties of grains in proto- 
planetary disks may provide important constraints on the 
formation history of circumstellar dust. 

We first focus on shape properties. Very little is 
known about the structure of grains in circumstellar disks. 
Particles may be compact or "fluffy" and may be chemi- 
cally homogeneous or consist of a mix of different mate- 
rials. Information may be obtained from measurements of 
grains sublimating from solar system comets nearing per- 
ihelion. From a comparison of measurements of the an- 
gle dependence of the scattering al bedo of several bright 
comets with theoretical predictions lOehrz fc Nevl <|l992h 



found the coma grains to be consistent with "fluffy" ag- 
gregates of smaller compact particles. Properties of the 
aggregate, such as temperature and spectroscopic signa- 
ture, will be affected by the actual degree of fluffiness of 
the grain, often e xpressed in terms of the poros ity fac- 
tor of the particle (jLisse et alJll998t lHarkerlll999j) . Large 
porous grains may heat to much higher temperatures than 
compact grains of the same mass and may show spectro- 
scopic characteristics similar to the smaller compact parti- 
cles constituting the fluffy aggregate. Theoretical consider- 
ations suggest cometary grains to be e xtremely fluffy with 
poros ity factors 0.93 ^ P ^ 0.975 l|Greenberg fc Haeel 
Il990h . effectively implying that some spectroscopically de- 
duced properties relate to the sm aller compact units m ak- 
ing up the porous aggregate fe.g. lBrucato et al"l fl999. and 
references there in) . In our modelling we will therefore con- 
centrate on (small) compact particles, leaving open the 
possibility that they may be part of larger fluffy aggre- 
gates. 

Re garding the chem i cal ho mogeneity of the small par- 
ticles, iLi fc Greenberd l(l998h pointed out that in situ 
mass spectra of comet P/Halley 1986 III dust showed 
that carbonaceous and silicate materials were mixed on 
fine scales, suggesting these two species are not physi- 
cally separated. If this is the case, the components are 
likely in thermal contact. The question of thermal con- 
tact is especially relevant with respect to the magnesium 
over iron content of the silicate material, a major con- 
stituent of circumstellar dust, as the absorption proper- 
ties of olivine (Mg2zFe2-2:cSi04, with x between and 1) 
and pyroxene (Mg^Fei^^SiOa) sensitively depend on this 
ratio. The pure magnesium silicates forsterite (Mg2Si04) 
and enstatite (MgSiOa) are optically much more trans- 
parent than iron rich silicates, implying that if Fe-rich 
and Fe-poor materials coexist as separate particles - with 
comparable shape and size - the Fe-rich dust will reach 
significantly higher temperatures. 

We have opted to treat the different chemical species 
as physically separated. In the case of Hale-Bopp this al- 
lows to investigate whether the forsterite crystals, of which 
the spectroscopic signatures are prominently visible in the 
mid-IR spectrum, are in thermal contact with the bulk 
amorphous silicate material, which is likely to contain 
a significant fraction of iron, by determining the mass- 
averaged temperatures of both species. If significantly dif- 
ferent, one may assume the components to be physically 
separated; if similar temperatures are found, it is likely 
they are in thermal contact. 

For each particle, we calculate the extinction prop- 
erties from optical constants determined in laboratory 
experiments, as listed in Table We assume spherical 
grains, for which we use Mie calculations, or a continuous 
distribution of ellipsoidal grains, for which we use CDE 
calculations, to determine the absorption and scattering 
coefficients (see iBohren fc Huffmanl l)l983|) for al full re- 
view on these methods). 

As we treat the CDE particles in the Rayleigh limit, 
we can not study grain size effects. For spherical particles 
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this assumption is not required. The spherical particles 
of our multi-component mixture of grains range between 
minimum size a m j n and maximum size a max and are dis- 
tributed following a power-law, i.e. 



n(o) oc 



(1) 



Theoretical calculations predict this type of size distribu- 
tion whenever there is shattering and coagulation of grains 
through grain-grain collisions (Bicrmann & Hartwit 1980). 
To get some feeling for the value of the power-law index, 
extinction observations imp ly a size distributio n with m 
= 3.5 for interstellar grains ijMathis et alJll977|h 



3.2. Chemical composition of grains 

The grain species used to model the spectra of both ob- 
jects are listed i n Table El The dust com position is similar 
to that used bv iBouwman et, al to model the iso- 

lated Herbig Ae stars AB Aurigae and HD 163296, though 
more appropriate sources for the optical constants of crys- 
talline silicates and water ice are used. Of the two iso- 
lated Herbig stars mentioned only in HD 163296 a small 
amount of forsterite could be identified, modeling of which 
did not sensitively depend on the adopted optical con- 
stants. However, for the two objects investigated here - 
which show prominent features of olivine - differences in 
laboratory measurements are important. To model the 
crystalline silicat e featu res we tried three sets of data: 
IServoin fe Pirioul l)l973|) meas ured th e optic al properties 
of forsterite (Mg 2 Si0 4 ; FolOOh lSteverl l)l974l) used a natu- 
ral olivine sample, which was estimated to contain a small 
amount of iron, i.e. x = 0.91 (Fo91). The sample used by 
iMukai fc Koike! (^990) contains an almost equal amount of 
iron (x — 0.90; Fo90). These measurements only cover the 
IR wavelengths. To estimate the optical properties of the 
grains at visual and UV wavelengths, where most of the 
stellar light is absorbed, we used measurements of compa- 
rable materials from several other sources, referenced in 
Table EJ 

Neither in AB Aur nor in HD 163296 the 44 /im feature 
characteristic for crystalline H2O ice was found, prompt- 
ing the use of amorphous water ice in modeling their 
circumstellar environment. In both HD 100546 and Hale 
Bopp, however, the 44 /iin feature appears to be visible, 
making it more appropriate to assume the water ice in 
these sources to be crystalline. 

4. Results 

4.1. HD 100546 

Plotted in Fig. El is our best model fit to the SED of 
HD 100546. The resulting model parameters are listed in 
Table The top panel shows the entire SED. A dashed 
line represents the Kurucz model for the central star. The 
lower panel shows the ISO-SWS and LWS wavelength re- 
gion. Indicated are the contributions to the spectrum of 



Table 2. Overview of dust species used. For each com- 
ponent we specify its latice structure, the wavelength in- 
terval over which optical constants are measured, and the 
destruction temperature. 



Species 


Latice 


Wavelength 


Tdcstr 


Ref. 




Structure 


H 


[K] 




[Mg,Fe]Si0 4 


A 


0.2-500 


1100 


(1) 


FeO 


C 


0.2-500 


1000 


(2) 


C 


A 


0.1-800 


1000 


(J) 


H2O 


C (ice) 


0.05-10 3 


150 


(4) 


Forsterite 


C 


0.04-3 


1400 


(5) 






3-250 




(6) 


Olivine (Fo91) 


C 


0.04-3 


1400 


(5) 






3-250 




(7) 


Olivine (Fo90) 


c 


0.01-0.3 


1400 


(8) 






0.3-2 




(9,10) 






7-200 




(11) 


Fe 


M 


0.1-10 5 


1500 


(12) 



Abbreviations used to designate the latice structure: A = 
Amorphous; C = Crystalline; M = Metallic. 

Refere nces: (1) |DorechjieX-et_al] lll995[): (2) [ Hen mng et "all 
ll995Tl: C3llPreibisch et alJ (Il993l): ?4llBertie et all lll969D: (5) 



IScottfe Dulevl Jl996l): r6nServoin fe Pirioul 1197S): (7) ISteveij 
Stapd Jl973|): (9) IJones fe Merrill! 

(12) lllenning fe'sTotrnienkoM ^ 



1 1974); (8) iHuffman fe Stapd il973T): (9) IJones fe Merrill 
1 1976ft: (10) iRoeers et alJ lll983ft fll) IMukai fe Koikd (Il990 " 



the individual dust components as listed in Table |3J To 
fit the SED a bi-modal grain size distribution is required, 
simila r as to that found in previous analysis of HAEBE 
stars l)Bouwman et al.ll200ot lBouwmanll200l|) . The small 
(< 10 /im) grains dominate the SED shortward of ~ 40/xm, 
while the large (up to 200 /im) grains dominate at the 
longest wavelengths. 

4.1.1. The small grain component 

Though the small grain component contains only a minor 
fraction of the total dust mass (~ 1%) it dominates the 
SED at near- and mid-IR wavelengths. To reproduce the 
observed fluxes in this wavelength range, a bi-model den- 
sity structure is required. As can be seen from Table 
the small grains are distributed between 0.3-9.8 AU, and 
between 9.8-43 AU, from here on referred to as zone 1 
and zone 2, respectively. The contributions to the SED of 
zone 1 and 2 are plotted in Fig. El The inner boundary 
of zone 1 at 0.3 AU is determined by the grain destruc- 
tion temperatures of the individual dust species as listed 
in Table El A similar dust component as that in zone 1 
is also present in ABAur, HD 163296 and HD 104247 
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Component: Small grains Large grains 



two component density structure 
zone 1 zone 2 



Mdust/Mg 




5.1 


10" 9 






6.7- 10~ 7 






6.5 


10~ 5 


Dust species 


R [AU] 


Pa [g 


r cm" 3 ] 


M frac 


R [AU] 


pa [gr cm -3 ] 


M frac 


a [/im] 


M frac 


a [^m] 


Amorph. silicate 


0.3-9.8 


1.6 


• lO" 17 


0.88 


9.8-43 


2.7- 10 -18 


0.7 


0.01-10 


0.53 


10-200 


Carbon 


0.3-9.8 


1.7 




0.09 


9.8-43 


3.1 • 10 -19 


0.08 


0.01-10 


0.05 


10-200 


Metallic iron 


0.3-4.9 


4.2 


■ io- 19 


0.01 








0.01-0.1 






Forsterite 


0.3-24 


6.0 


■ IO" 21 


0.02 


24-600 


1.4 • 10~ 21 


0.19 








Water ice 










9.8-43 


1.15 ■ 10~ 19 


0.03 


0.01-10 


0.42 


10-25 



Table 3. Best fit model parameters of HD 100546. Listed are the parameters defining the density and grain size 
distribution, the chemical composition and the mass fraction Mf lac of the individual dust species. Both the small grain 
and large grain component have a power law density distribution p{r) oc and grain size distribution n{a) oc a~ 2 . 
The radial extent of the small grain component is given in the table, large grains are present between 28-380 AU with 
a density at the inner boundary of po — 1-5 • 10~ 18 . Note that the densities given at the inner radius po — p(R; n ) 
of zone 1 assume all dust species to be present. However for temperatures greater than the destruction temperature 
Tdcstr (see Table EJ) the actual densities are slightly lower. 



l|Bouwman et al.ll2000l [Bouwmanll200l|) . The sudden in- 
crease in density at 9.8 AU, marking the onset of zone 2, is 
found only in HD 100546. It is this additional component 
that produces the much larger mid-IR luminosity com- 
pared to the other HAEBE systems as discussed in Sect.|21 
As can be seen from Fig. [3J the emission from zone 2 com- 
pletely dominates the SED between ~20 to 50 /on. The 
resulting total mass of 6.8 • 10~ 7 Mq in small grains is 
two to three orders of magnitude larger than found in the 
studies mentioned above. 

A major fraction of the dust emission seen in HAEBE 
systems is due to amorphous silicates. Though it is diffi- 
cult to determine the exact nature of the amorphous ma- 
terial, an excellent fit can be made if we use the optical 
properties of a silicate glass with an olivine stoichiome- 
try (see also Table 01. The grain size of the amorphous 
silicate is well constrained by the shape and strength of 
the 10 /im silicate feature, and the flux ratio between the 
10 and 18 /im bands. Apart from the amorphous silicates 
an additional source of continuum emission is required. 
We added this in the form of carbonaceous grains. Since 
this material has no distinct spectral features, the grain 
size and consequently, mass fraction, are less well estab- 
lished and are uncertain within a factor of two. In our 
best fit model we assumed an equal grain size distribu- 
tion as for the amorphous silicate. The band at 43.3 /tm 
in the spectrum of HD 100546 indicates the presence of 
a population of small ice grains. Given the ice sublima- 
tion/evaporation temperature of 150 K, water ice can only 
be present in zone 2 outside of the "snow limit". For the 
small ice grains we also assumed the same grain size dis- 
tribution as for the amorphou s material. Interesting to 
note is that ISmith et al 1 lll994h showed that the ice band 
can only shift to the observed wavelength of 43.3 /jm if 
amorphous ice grains where heated to the crystallization 



temperature and then cooled to ^50 K. One could inter- 
pret this as a local (i.e. within the "snow" regime) heating 
event or local radial mixing. 

As one can see from the inset in Fig. |2J the near-IR 
flux shortward of 4 /im is dominated by the emission from 
metallic iron grains. These grains are the only refractory 
dust species which are stable at the required high tem- 
peratures (~ 1500 K), and that have sufficient opacity at 
near-IR wavelengths to emit in this spectral region. We 
had to constrain the spatial distribution of the metallic 
iron grains to within 4.9 AU. The grain temperature of the 
iron grains is between 1500 and 750 K. This suggests that 
metallic iron is formed at high temperatures, most likely 
by solid state reduction of iron bearing silicates heated 
above the glass temperature at the very inner parts of the 
circumstellar disk. 



4.1.2. The crystalline silicate component 

In the sample of isolated HAEBE stars presented in 
iMeeus et"alTl)200l|) . HD 100546 shows the most prominent 
features of crystalline silicates. This enables one, at least 
in principle, to determine the exact chemical composition 
and crystalline structure of the silicates. We used three 
sets of laboratory data listed in Table to determine the 
best fit to the crystalline silicate bands. Plotted in Fig. [31 
are the 24 and 34 /xm silicate bands of HD 100546, com- 
pared with our best fit models for the three data sets. As 
o ne can clearly see, the mo del based on the measurements 
of lServoin Pirioul l|l973h provides the best match with 
the ob s ervati ons. For this model and the model using the 
ISteverl l|l974h data set, the CDE shape distribution (see 
Sect. |3J) gives a much better fit than do sphe rical grains. 
In the model using the lMukai fc Koikel l|l990h data, MIE 
calculations produce the best match. We conclude that 



Bouwman et al.: HD100546 



7 




10 100 

Fig. 2. The top panel shows our best model fit to the spec- 
tral energy distribution of HD 100546 (thick solid line). 
Indicated with the dashed line is the spectrum of the cen- 
tral star. Triangles indicate ground-based and IRAS pho- 
tometry. Also plotted are the the ISO-SWS and LWS spec- 
tra. The dashed line, the thin solid line and the dashed- 
dotted line represent the contributions of the small grain 
components in zone 1, zone 2 and the large grain compo- 
nent, respectively. The bottom panel shows our fit to the 
ISO spectra in detail. Also plotted in this figure are the 
relative contributions of the individual dust components. 
Indicated with a solid line are the contributions of amor- 
phous olivine and forsterite marked in the figure with 01 
and Fo, respectively. The contribution of water ice, car- 
bonaceous grains and metallic iron are represented by the 
dashed, dotted and dashed-dotted lines, respectively. The 
curves of the carbonaceous, forsterite and metallic iron 
grains are offset by -40, -75 and -90 Jy, respectively. The 
inset shows the 2 to 8 fim region. 

the crystalline silicate features seen in the ISO spectra of 
HD 100546 are due to non-spherical forsterite grains. In 
our best fit, presented in Fig. [5] and Tabled we used this 
latter model for the crystalline silicate. 

The forsterite grains also require a bi-modal density 
structure (zone 1 and zone 2). As can be seen from Tabled 
the crystalline mass fraction is much higher in zone 2 than 
in zone 1, indicating that the abundance of crystalline ma- 
terial is increasing outwards. In our best fit model the 
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Fig. 3. The 24 and 34 fim forsterite bands of HD 100546. 
O yerplotted are the best fit silicate b ands us i ng the data 
of lServoin fc Pirioul lll973l) (solid lin e'l. lsteverl lll974h (dot- 
ted line) and lMukai k. Koikd l)l990|) (dashed line). In the 
latter model spherical grains where assumed with a size 
of 0.01 /im. The other models have a CDE shape distribu- 
tion. 



forsterite grains are not co-spatial with the other dust 
species. The spatial distribution of forsterite grains has to 
extent out t o 600 AU, well beyond the observed disk limit 
of 380 AU l|Augereau et all 1200 1|) . This, however, could 
be an artifact of the grain model we used for calculating 
the optical properties of the forsterite grains. The CDE 
approximation assumes the grains to be in the Rayleigh 
limit, and consequently the optical properties of the par- 
ticles do not depend on grain size. However, outside the 
Rayleigh limit, particle size has a strong effect on grain 
temperature. Larger grains have to be placed much closer 
to the central star than smaller grains to reach the same 
equilibrium temperatures. The grains can however not be 
much larger than 1 /im to be able to reproduce the ob- 
served spectral features. This is smaller than the typical 
grain size of the amorphous material. Our best estimate 
for the effect of grain size on the spatial dist ribution, using 
the measurements of lMukai fc Koikd l)l990j) indicates that 
the crystalline grains could be moved inwards to within 
the observed disk limits, but not far enough to be co- 
spatial with the amorphous silicate. 

The grain size distribution also influences the mass 
fraction of the dust species. The model assuming a max- 
imum grain size of 2 /xm for the forsterite grains would 
increase the mass in crystalline dust by a factor of five. 
This would imply that as much as 50 % of the total mass 
of the small grain component could be crystalline silicates. 
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4.1.3. The large grain component 

To reproduce the observed slope and fluxes at mm wave- 
lengths a dust component with a large grain size up to 
200 /im is required. Their contribution to the SED of 
HD 100546 can be seen in Fig. Due to this large gram 
size, one can not determine its exact composition from ISO 
spectroscopy, since such large grains will radiate as black 
bodies at these wavelengths. We used the dust species 
found spectroscopically to model the large grain compo- 
nent and assumed their relative mass ratios to be identi- 
cal to that of the small grain component. The mass frac- 
tions of the large amorphous silicate and carbonaceous 
grains, listed in Tabled are therefore uncertain. Similarly, 
it is not possible to determine if these large grains contain 
crystalline material as this will also result in black body 
emission. We therefore did not use crystalline silicates in 
the large grain component. However, due to the presence 
of the emission band at 60 /im attributed to crystalline 
water ice grains, the mass fraction of this dust species 
can be determined. The maximum size the ice grains can 
have and still produce a 60 /im band, is 25 /im. This is 
much smaller than the maximum size required to fit the 
slope of the SED at mm wavelengths. There are however 
some problems with the 60 /im band as we will discuss in 
Sect. gHgl 

We have constrained the spatial ext ent of the large 
grains using the disk size observed by lAugereau et alJ 
l|200l|) with the HST/NICMOS2 instrument. To fit the 
observations, a gap of 28 AU is required between the cen- 
tral star and the dust shell containing the large g rains. 
As noted in previous studies l|Bouwman et al-lfeoOOj) . this 
gap is most likely an artifact of the assumption made in 
this study that the medium is optically thin in the radial 
direction. In HD 1005 46 we see the disk fairly face on (51° 
lAugereau et al.ll200lj) . implying that in the line-of-sight 
the medium may well be optically thin. Therefore, the de- 
rived mass over temperature distribution T^m) remains 
m eaningful (albeit not th e ?d(r) structure). As suggested 
bv lBouwman et aTl |2000) the large grains have most likely 
settled to the disk mid-plane, extending all the way to the 
disk in ner edge. Recent full 2D radiative transfer calcu- 
lations (|Bouwmanll200lj) confirm this qualitative picture. 
These calculations also show that the ?d(r) of the hot 
small grains, residing in the optically thin surface layer of 
the disk, is in agreement with the simple spherical opti- 
cally thin model adopted in this paper. 



4.1.4. The mass temperature distribution 

The SED as presented in Fig. |3 is determined by the 
mass-temperature distribution of the circumstellar dust. 
Irrespective of the assumed model geometry, this distri- 
bution has to be reproduced for dust that is contributing 
to the observed emission features. Plotted in Fig. [5] is the 
derived mass over temperature distribution of the best fit 
model. The upper two panels show the mass-temperature 
distribution of the small grain component, the lower panel 




X [|lm] 

Fig. 4. Contributions to the model SED of the individual 
dust species per zone. The first two panels show the contri- 
butions of the small grain component with sizes < 10 /mi 
located in zone 1 and zone 2, respectively. The lower panel 
shows the distribution of the large grains (> 10 /im). The 
lines represent the same materials as used in the lower 
panel of Fig. [21 The top solid lines indicate the contribu- 
tions of amorphous silicates and the lower solid lines the 
contribution of forsterite. 



shows the same for the large grains. Indicated in the figure 
are the individual contributions of each species. The ver- 
tical axis unit is chosen in such a way that the integral of 
the cumulative mass-temperature distribution equals the 
total dust mass. We have plotted in Fig. 0] the relative 
contribution to the SED of the individual dust species to 
get a feeling of their importants to the overall fit. Plotted 
in Fig. El for comparison is the mass-temperature distri- 
bution of ABAur as derived bv iBouwman et al] |2000). 
The first panel shows the mass-temperature distribution 
of the small grains, the second panel that of the large grain 
component. 

A comparison between Figs. and HO shows that the 
dust in HD 100546 is characterized by a bi-modal mass- 
temperature distribution, similar to that in the AB Aur 
system. The inset in the top panel of Fig. [5] shows the 
mass distribution at temperatures above 500 K, produc- 
ing the near-IR emission. This component, dominated by 
hot iron or iron oxide grains, is very similar for both ob- 
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Fig. 5. Cumulative dust mass 
over temperature distribution of 
HD 100546. Indicated in the panels 
are the relative contributions of the 
individual dust species. The first 
two panels show the contributions 
of the small grain component with 
sizes < 10 /zm located in zone 1 
and zone 2, respectively. The lower 
panel shows the distribution of the 
large grains (> 10 The inset 

in the top panel shows the high 
temperature dust component of 
HD 100546. 



jects. iMeeus et al.l l|200l|) have shown that the near-IR 
spectrum of HAEBE systems is generic, suggesting a sim- 
ilar geometry and/or dust composition at the very inner 
parts of the circumstellar dust region. However, the total 
mass contained in the small grains in HD 100546 is almost 
three orders of magnitude larger. This excess mass is lo- 
cated in zone 2 (see Table|3J) of the small grain component, 
which is not present in AB Aur and is the cause for the 
much larger mid-IR excess as observed in HD 100546. 

Focusing on the mass-temperature distribution of the 
small grains, one notices the difference in the distribution 
between the amorphous dust and the forsterite grains. To 
fit the SED, the crystalline silicate must have a much lower 
temperature than the amorphous grains. This result holds 
irrespective of the adopted laboratory measurements of 
the crystalline silicate. This mass-temperature distribu- 
tion imposes the small forsterite grains to have a larger 
radial extent compared to the amorphous dust. 



The large grain component of HD 100546 has a distri- 
bution similar to that of AB Aur. The total mass contained 
in this component is the same within a factor of two. The 
small difference in dust mass between both systems is due 
to a grain size difference, where in AB Aur the maximum 
grain size is smaller by a factor of two. 

4.1.5. The residual spectrum 

Plotted in Fig. is the residual spectrum of HD 100546, 
after subtraction of the best fit model to the ISO-SWS and 
LWS spectra. Shown in the top panel is the 2 to 15 /jm re- 
gion in which the emission bands from polycyclic aromatic 
hydrocarbons (PAHs; marked with the vertical lines) are 
clearly visible. For comparison we also plotted the PAH 
spectrum of the [WC] star BD+30 3639, which may serve 
as a template of a generic PAH spectrum. Apart from 
the bands at 6.25, 7.9 and 11.3 ^m, also PAH features at 
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Fig. 6. Cumulative dust mass over 
temperature distribution of AB Aur 
and comet Hale-Bopp. The two top 
panels show the distribution for 
AB Aur for small (< 10 /im) and 
large (> 10 /im) grains, respectively. 
The bottom panel shows the distri- 
bution for the solar system comet 
Hale-Bopp. Indicated in the panels 
are the relative contributions of the 
individual dust species (see Fig. El 
for identification). 



longe r wavelengths can be observed (see also iHonv et alJ 
l200llL 

The lower panel shows the 35 to 150 /im region. 
Marked with shaded areas are residual emission bands. 
The band at ^100 /im has bee n interpreted as evid ence for 
the presence of hydro silicates l(Malfait et alJll99^ . If cor- 
rect, this can be seen as evidence for the presence of larger 
bodies on which the hydro silicates are believed to have 
formed by aqueous alteration of liquid water. For compar- 
ison we plotted the solid-state band of HD 142527 which is 
identified with hydro silicates (Mal fait et al.ll999l) . As one 
can see the width and position are similar for both bands. 
However, a secure identification of this rather broad fea- 
ture strongly depends on the exact run of the underlying 
continuum. A small change in the location of this con- 
tinuum may change both the strength and width of the 
feature or may even cause it to disappear. 



The residual at ~55 /im could indicate the presence of 
an emission band of a species not treated in our models. It 
is, however, more likely that this residual originates from 
the poor fits of the ice bands at 43 and 60 /im. Though 
the positions of the observed ice features coincide well with 
our ice model, the width and relative strength of the 43 
and 60 /im bands do not. This could reflect limitations of 
our simple ice model, e.g. the assumptions of homogeneous 
spherical grains. Also the temperature of the material at 
which the optical properties where measured (100 K for 
our data) will influence the strength and width of the fea- 
tures. The observed position of t he ice band at 43. 3 /im 
indicates that the ice is ^50 K l|Smith et al.l ll994K An 
additional problem is that the ice features appear in two 
separate instruments, LWS and SWS. These instruments 
connect at 45 /im (indicated with an arrow in Fig.|7|). The 
error in the absolute flux calibration between both instru- 
ments is 10-15 %. Lowering the LWS spectrum with this 
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Fig. 7. Residual of the ISO-SWS and LWS spectrum of 
HD 100546 after subtraction of our best fit model. The 
top panel shows the 2 to 15 /im region clearly showing 
the PAH emission bands (indicated with vertical lines). 
For comparison we also plotted the PAH spectrum of 
BD+30 3639 (offset by -12 Jy). The bottom panel shows 
the 35 to 150 fim region. Indicated with an arrow is the 
connection between the ISO-SWS and LWS spectra. Also 
plotted in the figure are the absorption spectra of several 
silicates of interest multiplied by a 70 K black body, and 
the emission band in th e spectrum of HP 1 42527 identi- 
fied with hydro silicates {Eylalfait ct al. 1999). The shaded 
areas indicate features in the residual spectrum, discussed 
in Sect. l4~l~5l 



percentage would enable us to fit the ice bands without 
the residue at ~55 ^m. Interesting to note is that the LWS 
spectrum would then also agree better with the IRAS pho- 
tometry which is now below the LWS spectrum, as can be 
seen from Fig.El However, we stress that from the ISO ob- 
servations there is no strong evidence for a jump between 
the ISO-SWS and LWS spectra at ~45 /mi. 

A third feature in our residual spectrum appears at 
40.3 /im. This band could be due to crystalline dust species 
other than forsterite. Plotted in Fig. [7| are the spectra 
of di opside and ortho- and clino-enstatite l|Koike et alJ 
200CJ) . All these species show bands in the 40 /im re- 
gion. The 40 /im band of diopside coincides nicely with 
the residual spectrum but would also produce a band at 
65 /xm at the location of the water ice band, making the 
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Fig. 8. Best model fit to the ISO-SWS spectrum of Hale- 
Bopp. Also plotted in the figure are the contributions of 
the individual dust species to the spectrum. The solid lines 
show the contribution of the amorphous silicate with an 
olivine stoichiomctry (marked in the figure with 01) and 
forsterite (Fo). The dotted, dashed and dot-dashed lines 
show the contributions of carbonaceous material, ice and 
iron oxide grains, respectively. The curves of forsterite, 
carbonaceous, ice and iron oxide are offset by -10,-50,-70 
and -90 Jy, respectively. 



problem at 55 /im even more severe. If one of the suggested 
materials is responsible for the 40 /im band, it has to be 
cold, i.e. Si 70 K, as at shorter wavelengths no evidence 
for crystalline silicates other than forsterite can be found. 
Finally we mention that the bands at 40.3 and 43.3 /im, 
coinciding with the ice band, may be due to forsterite. 
This suggestion is bas ed on mass absorption spectra taken 
bv I Jaeer et al J 1^998), which show evidence for bands at 
these wavelengths. Though these features are weak, given 
the overall strength of the forsterite bands this could be a 
possibility. 



4.2. Hale-Bopp 

Plotted in Fig.|5|is the ISO-SWS spectrum of comet Hale- 
Bopp, together with our best fit model. Also shown are 
the contributions to the spectrum of the individual dust 
species. The parameters of the model fit are listed in 
Table IH As has be en discussed by previous authors (e.g. 
iMalfait et alll998|) . the ISO spectrum of Hale-Bopp bears 
a striking resemblance to that of HD 100546. Our modeling 
efforts of the spectrum of Hale-Bopp will focus on estab- 
lishing the similarities and differences in dust properties 
between the comet and HD 100546. 
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Table 4. Model fit parameters of Hale-Bopp. Listed are 
the parameters defining the density and grain size dis- 
tribution, the chemical composition and the mass frac- 
tion Mf rac of the individual dust species. For all dust 
components we assumed a power law density distribution 
p(r) oc r^ 1 and size distribution n(a) oc a -2 8 . The radial 
extend and grain size range are given in the table. 



4.2.1. Grain composition and properties 

In general we find that the dust composition of Hale-Bopp 
is similar to HD 100546, with a few exceptions. We find no 
evidence in the spectrum of Hale-Bopp for emission from 
metallic iron grains as in the case of HAEBE systems. 
This could be a spurious result as in the cold (160 K) 
dust of Hale-Bopp iron grains have no distinct spectro- 
scopic feature and their emission could blend with that 
of other smooth continuum sources, for instance from car- 
bonaceous grains. However, it could also imply that metal- 
lic iron grains are really absent, given the suggested high 
temperature formation mechanism of this grain species 
(see Sect. |4~TTT|) . At lower temperatures meta llic iron ca n 
also react to form iron sulfides or oxides fe.g lGaillll998|) . 
Indeed, we find evidence for emission from these latter 
components, in this paper modeled with thermal emission 
from iron oxide grains, similar to tha t found in AB Aur 
and HD 163296 l)Bouwman et al.ll2000|L though less abun- 
dant. The possibility, however, exists that the emission 
component here attributed to iron oxide is (partially) due 
to iron sulfide whic h has similar emiss ion properties in this 
wavelength region l|Keller et al.ll2000ll . 

A component of amorphous silicates and carbonaceous 
grains similar as in HD 100546, constituting a major frac- 
tion of the dust mass, is seen in Hale-Bopp. Interesting 
to note is the typical grain size of the cometary dust, 
which follows a power-law slope m = 2.8, much larger 
than in HD 100546 (m = 2). Note that for the forsterite 
and iron oxide grains we used a CDE shape distribution 
which does not allow for the determination of grain size. 
Spectroscopically, the best fit to the forsterite bands is, 
li ke for HD 100546. achiev ed by using the measurements 
of lServoin k, Pirio u Hl973h . Similar uncertainties apply for 
the derived forsterite mass fraction as for HD 100546, dis- 
cussed in Sect. l4.l!2l Our best estimates for the effects of a 
grain size distribution on the derived mass indicates that 
the maximum mass fraction of forsterite could be about 



"40 %, i.e. six times higher than in HD 100546. A grain size 
distribution would also result in a higher continuum con- 
tribution of the forsterite grains, which eliminates the re- 
quirement for an additional component in the form of iron 
oxide. From the ISO-SWS spectrum it is difficult to deter- 
mine the characteristics of t he water ice componen t. We 
therefore used the results of iLellouch et alJ l|l998f) . who 
found from an analysis of the ISO-LWS spectrum that the 
average water ice gra in is 15 um. The mass fraction of the 
ice is similar to what ILellouch et alJ 1 1998h find, however, 
they infer a total dust mass of 1.1 • 10 11 kg, almost two 
orders of magnitude larger than we derive. They however, 
assumed a grain size of 100 /im in deriving the total dust 
mass. Lowering this to the sizes we find would scale their 
total dust mass to a similar value as found here. 

Previous studies (e.g. IWooden et al.l fl999^) have also 
used pyroxenes to model the observations of Hale-Bopp. 
However, from the ISO spectrum no evidence can be found 
for the precence of such a component. We have therefore 
not used these dust species in our modeling. 



4.3. Thermal contact of grain species 

The lower panel of Fig.EJshows the mass-temperature dis- 
tribution of Hale-Bopp. Contrary to the HD 100546 sys- 
tem, the dust species in Hale-Bopp all have the same mass 
averaged temperature of approximately 160 K. The differ- 
ence in the width of the distribution of the individual dust 
species results from the assumption of a grain size distri- 
bution for the amorphous dust, giving a broader distribu- 
tion than for the other dust species (forsterite, iron oxide) 
where a CDE shape distribution is assumed. In order to 
model the spectrum of comet Hale-Bopp, we found, as ex- 
pected, that the bulk of the dust material (i.e. the amor- 
phous silicates and the carbon, together accounting for 87 
percent of the dust mass) has to be placed at the proper 
heliocentric distance at the time of observation (2.8 AU 
and distributed over 0.2 AU), giving it the required mass 
averaged temperature of ^160 K. To fit the near-IR wave- 
length range, positioning the iron oxide at the correct dis- 
tance resulted in a temperature for this species too high 
to properly fit the spectrum. We determined the best fit 
temperature for this material by positioning it at a larger 
distance. This yields T ~ 160 K. The same approach was 
followed to determine the temperature of the water ice and 
the crystalline silicate. Here it was found that due to their 
relatively poor light absorption properties both materials 
had to be placed closer to the star, but again at such a dis- 
tance that their bulk mass averaged temperature is 160 K. 
A co-spatial distribution of all grain species, complying 
with the derived mass-temperature distribution and the 
observational constrains on the heliocentric distance, can 
only be obtained if thermal contact between the grains 
is imposed. This result implies that all dust particles are 
in thermal contact with one another and that their tem- 
perature is essentially determined by the materials domi- 
nating the thermal energy content of the dust. The latter 
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are the amorphous silicates and the carbonaceous mate- 
rial. The typical dust temperature we find is much lower 
than t hat addopted in previous studies bv lWooden et all 
l)l999(l . They find a typical temperature for the amorphous 
silicates and crystalline olivine of 245 K, which is ~ 85 K 
higher than what we find. Also, they suggest a temper- 
ature differnce between crystalline pyroxene grains and 
other dust species which is contradictory to our findings. 

Though the dust species are the same as in HD 100546, 
Fig. El shows the fundamental difference between the dust 
in Hale-Bopp and that in the Herbig Be system: The 
grains in Hale-Bopp are in thermal contact, those in 
HD 100546 are not. Forcing the grain species in HD 100546 
to be in thermal contact, the observed difference in the 
mass-temperature distribution between the amorphous 
and crystalline dust can not be reproduced. This leads 
to the conclusion that the grains in HD 100546 are homo- 
geneous, separate entities. 

5. Discussion 

To summarize, the main results from the analysis of the 
SED of HD 100546 are the following: (1) The SED of 
HD 100546 stands out compared to other well studied iso- 
lated HAEBE systems in the sense that it has a much 
weaker near-IR excess, but a much stronger mid-IR ex- 
cess. (2) An additional component of small (<10 /urn) dust 
grains is present in the circumstellar material surrounding 
HD 100546, not seen in other systems. (3) HD 100546 has 
by far the largest crystalline silicate mass fraction of all 
HAEBE stars studied with the ISO satellite. (4) A com- 
parison of the mass-temperature distribution of the small, 
spectroscopically identifiable amorphous and crystalline 
silicates shows that the bulk of the crystalline silicates 
has a much lower temperature than that of the amorphous 
material. From these results several conclusion concerning 
the disk geometry and the origin of the crystalline silicates 
can be drawn, as we will discuss below. 

5.1. Disk Geometry 

As we showed in Sect.0 the SED of HD 100546 deviates 
considerably from that of ABAur. The SED of the lat- 
ter can be reproduc ed with a hydrostatic r adiative equi- 
librium disk model IjDullemond et al.ll200l|) . The unusual 
luminosity distribution of the circumstellar material as 
listed in Table ^ however, can not be reproduced with 
a model having a similar disk geometry as used to model 
ABAur. For dust grains to efficiently re-radiate the ab- 
sorbed stellar light in the mid-IR, the typical dust tem- 
perature must be in the order of 200 K, which is the equi- 
librium temperature of grains in the disk surface layer 
at a radial distance of about 10 AU. To reproduce the 
observed mid-IR luminosity of 0.35 the relative scale 
height H/R of the disk surface at the radial distance of 
~10 AU has to be 0.35 in order to "capture" this amount 
of starlight. This is a factor of two higher than the de- 



rived scale height at this radius of the d isk surface in the 
ABAur system l|Dullemond et al.ll200l|) . 

However, as iDullemond et al ] ll200lh show, if a hole 
forms in the disk, the inner edge of the disk will "puff up" 
as a consequence of the direct irradiation by the central 
star. The scale height at the edge will typically increase 
by a factor of two over the height of a normal flaring disk 
without inner hole. Using the same model parameters as 
derived for AB Aur, but with a hole of 10 AU would re- 
produce the required disk scale height and consequently 
the observed mid-IR excess of HD 100546. However, the 
region inside 10 AU cannot be free of dust since there is 
still excess emission in the near infrared. The material pro- 
ducing this excess must be located at distances of ^1 AU. 
Instead of an inner hole, the disk should therefore have 
a gap starting outside of ~ 1 AU and extending to about 
10 AU. The question is then, whether enough stellar radi- 
ation can reach the inner boundary of the outer disk in 
order to still cause the i ncreased scale height there. In the 
Dullemond et alJ (|200lJ) model, the material in the inner 
regions casts a shadow over parts of the disk. However, as 
we have shown in Sect. 2, the near-IR excess of HD100546 
caused by the inner part of the disk is relatively small, a 
disk surface height of H/R = 0.085 is geometrically suffi- 
cient to absorb the required 8.5% of the stellar radiation. 
Therefore the blocking of light by the hot dust is only a 
small effect in HD100546. We can now identify the mid-IR 
component in the ISO spectra with a "wall" , being the far 
side of a disk gap directly illuminated by the star. The 
additional dust component of small grains (zone 2, see 
Table EJl, unique to HD 100546, should then be the dust 
being present in this "wall" at ~10 AU. 

5.2. Minearology: the origin of the crystalline silicates 

As was already noted in previous studies, HD 100546 
has the highest k nown fraction of crystalline silicates of 
all HAEBE stars llMalfait et al.lll99St iMeeus et al1l200lt 
iBouwman et alJ l200l[) . More puzzling, the bulk of the 
crystalline silicates is at low temperatures (~ 100 K, see 
Fig. EJ). The question is how to explain both the high 
mass fraction as well as the mass-temperature distribu- 
tion of the crystalline silicates. Fig. [5] shows the relative 
forsterite mass fraction, as plotted in Fig. [5] relative to the 
total mass in silicates as a function of temperature. Clearly 
visible is the rise in the crystalline mass fraction at lower 
temperatures. A crucial question is which mechanism is re- 
sponsible for this outward increase of the crystalline dust 
fraction. This question is closely related to the formation 
processes of the crystalline silicate. Forsterite will form 
by annealing of amorphous silicate at temperatures above 
the glass temperature. At temperatures above 1100 K the 
annealing time-scale b ecome much smaller th an any dy- 
namical time-scale fe.g. lHallenbeck et al.l200o|) . and a fast 
transformation into crystalline silicates follows. This im- 
plies that the dust in the very inner parts of the circum- 
stellar disk, with temperatures above 1100 K, will be en- 



14 



Bouwman et al.: HD100546 



1.00 



g 

"o 
ctf 



co 
co 
re! 

E 

CD 



CD 
-i— > 

co 
o 



0.10 




0.01 L. 
3.5 



3.0 



2.5 
logT 



2.0 



1.5 



Fig. 9. Relative forsterite mass fraction as a function of 
temperature. The solid line indicates our best model fit. 
The dotted line represents a model where the crystalline 
and amorphous silicates are co-spatial with a constant rel- 
ative mass fraction as a function of radial distance. Note 
that between the mass-temperature and spatial density 
distributions an opacity effect exists. The forsterite grains, 
having lower opacities in the UV and optical compared to 
the amorphous silicate, have a different temperature dis- 
tribution. This shift in the temperature distribution be- 
tween the amorphous and crystalline dust causes the rise 
seen around ~300 K for the model with constant relative 
mass fraction as a function of radial distance. 

tirely crystalline up to the evaporation temperature of sil- 
icate grains at about 1400 K. Several radial mixing mod- 
els have been proposed to transport material outwards 
to the regime where crystalli ne silicates are ob served well 
below the glass temperature. IStevensonl l|l99fl() predicted 
that the mass fraction of the mi xed ma t erial w ould decline 
outwards following a power law. IPrinnl ((1990) argued that 
the neglect of nonlinear momentum terms in the Stevenson 
model could underestimate the radial mixing. He predicts 
that the mass fraction of the mixed material could almost 
be constant as a function of radius. Similar resu l ts are 
found in recent models bv lBockelee-Morvan et"all (|2000j) 
and lGaill (|200l[) . Plotted in Fig. ED (dotted line) is a model 
with a constant relative mass fraction of forsterite over 
amorphous silicate as a function of radial distance, repre- 
senting the case of the most effective radial mixing pre- 
dicted. As can be seen from this figure, outward radial 
mixing of material from the inner parts of the proto- 
planetary disk can not be the explanation for the observed 
distribut ion of crystalline silicat es. 

Also. lBouwman et al.l l)200lh show that the dust com- 
position producing the 10/im silicate feature, probing 



grains at the inner parts of the proto-planetary disk, is 
inconsistent with thermal annealing. Clearly a different 
mechanism for the formation of the observed crystalline 
dust fraction and distribution must be formulated. Having 
ruled out radial mixing, this implies that the crystalline 
silicates must have formed locally i.e. further out in the 
disk where they are presently observed. What could have 
caused this crystallization in the cool outer regions of the 
disk? One way of providing enough energy for the silicates 
to crystallize is to incorporate them into larger bodies. 
Radioactive decay of species like 26 Al could cause the (par- 
tial) melting and consequently crystall ization of silicate s 
in sufficiently large parent bodies (e.g. iHuss et al-llioOlj) . 
Afterwards, small crystalline dust grains can be produced 
by collisions between the processed parent bodies. If the 
parent body was large enough to fully melt during its ac- 
cretional phase, it will be differentiated. The crust will be 
deprived of siderophile elements like Fe/Co/Ni. The com- 
position of the crust will therefore largely be a magnesium- 
rich silicate, very similar to what is indicated by observa- 
tions. For a collisional cascade of initially large objects 
to proceed, the orbits must be sufficiently non-circular to 
produce large relative velocities. For this one needs many 
Pluto sized objects or a large Jupiter type p lanet (e.g. 
ICharnoz et al lEfioTl iKenvon fc Bromlevll200lh. 

An other possible mechanism for in s itu production of 
crysta lline silicates has been proposed bv lHarker fe Deschl 
( 2002). They put forward the possibility that shock waves 
triggered by gravitational instabilities could heat small 
amorphous grains to sufficiently high temperatures to 
crystallize. However, this mechanism still implies thermal 
annealing of small dust grains, and the 8.6 /im feature of 
silica should still be produced in this case. 

5.3. Proposed scenario 

As we have discussed in the previous two subsections, 
there are two things that set HD 100546 apart from other 
HAEBE stars, being its disk geometry and its mineralogy. 
We propose a scenario which can explain and connect both 
deviations. In this section we will speculate that a massive 
planet has formed in the disk at a distance of about 10 AU 
from the star, and that this planet is responsible for the 
striking difference in both the mass-temperature distribu- 
tion and in the chemical composition and crystallographic 
structure of the dust around HD 100546. 

A massive planet will be able to strongly influence 
planetesimals and the structu re of the g as disk within its 
Hill sphere which is given by l)Hil]Hl878li 

where A is the semi-major axis of the planets orbit, M p \ 
is the mass of the planet and M* is the mass of the star. 
iLin fe Papaloizoul l|l993j) showed that a planet will open 
a gap in the disk if its Hill sphere equals or exceeds the 
vertical disk scale height. Assuming that before the for- 
mation of a planet the disk was very similar to that of 
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AB Aur, the relative disk scale height at 10 A U would have 
been H/R = 0.045 ijDullemond et alJl200l|) . Substituting 
this into Eq. @, taking into account that the far side 
of the gap is "puffed up" by a factor of about two (see 
Sect. |21 and Sect. 15. If) , would yield a planetary mass of 
Mpi > 5.6Mj, indicating a Jupiter type planet. A forma- 
tion distance of ~10 AU for a gas giant seems reasonable 
to assume, in particular since the temperatures in the disk 
around a late type B star will be similar at 10 AU to 
the conditions at 5 AU in the early solar system. It has 
been suggested that the snow-line, i.e. the line where wa- 
ter ice can condense and remain stable in the disk, has 
enhanced the densities in the solar disk locally to spee d 
up the accretion of Jupiter ijStevenson &: Luninelll988|) . 
This ice condensation zone for the HD 100546 system is 
located at about 10 AU. Also, the time scale for the for- 
mation of a giant planet in the disk regi on is a few Myrs 
^Pollack et al.lll996t IrTersant et al.ll200l|) . which can eas- 
ily be accou nted for given the estimated a ge of HD 100546 
of >10 Myr l(van den Ancker et alJll998j) . A planet open- 
ing a gap in the disk would also be a natural explanation 
for the lower near-IR excess of HD 100546 compared to 
AB Aur. Such a gap could cause the depletion of the inner 
disk by starv a tion o f material from the outer disk, though 
iLubow et al.l l|l999fl argue that for planets of ~Mj still a 
net flow through the gap can exist not completely clearing 
the inner region. 

As briefly mentioned in Sect. 15.21 the formation mech- 
anism for the crystalline silicates could be a collisional 
cascade of planetesimals. A Jupiter mass planet can easily 
"heat up" the orbital elements of a swarm of planetesimals 
through gravitational stirring, causing sufficiently high rel- 
ative velocities for collis ions to be destructive. Indeed, 
IStern fc Weissmanl l|200l|) show that local collisional ero- 
sion of planetesimals may be a more important process 
for planetesimals near a giant planet than ejection to the 
Oort cloud. An important question is now, how much re- 
processed dust can be expected from a collisional cascade 
of such bodies. iDohnanvil Jl969) has shown that the col- 
lisional cascade leads to a size distribution f(a) oc a~ 3 ' 5 
or /(to) oc to -11 / 6 in which the mass is dominated by the 
largest particles. If the observations measure all particles 
up to a certain size limit a\ ~ 10 /im and the maximum 
size of planetesimals involved in the collisional cascade is 
100 km, then one finds m v i s /mtot = 10~ 5 . For the mea- 
sured amount of crystalline dust of m v - ls = 1.3 x 10~ 7 M , 
one finds a total mass of m to t = 1-3 x 10~ 2 M Q . This is 
more than there is available in the disk. Therefore, the 
observed dust cannot be the steady state solution of a 
collisional cascade. 

However, the steady state solution assumes that small 
particles are constantly removed at the low end of the size 
distribution. If the disk around the star were already gas- 
free, particles smaller than about 10 urn in size would b e 
blown away by radiation pressure fe.g. lArtvmowiczll988|) . 
However, the disk around HD 100546 very likely is not 
gas free. Even if large bodies can move around freely on 
purely Keplerian orbits, the motions of small grains will 
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Fig. 10. Schematic representation of the HD 100546 sys- 
tem. A giant planet has formed a gap in the flaring disk 
at 10 AU, causing an increase of the disk scale height. 
A collisional cascade induced by the giant planet is pro- 
ducing dust and throwing it up out of the disk mid-plane. 
The small crystalline grains above the disk are blown away 
by radiation pressure and rain down further out at a few 
hundred AU on the surface of the disk (arrows). 



be slowed down. This will keep small grains in the disk, 
and it will also drastically reduce the collision velocities 
between small grains. Therefore, small grains are no longer 
destroyed in collisions. Grains which are small enough will 
accumulate in the disk. The density of these grains will 
not be in collisional equilibrium. It will instead mirror the 
entire mass of larger bodies destroyed earlier. This greatly 
softens the mass constraint we derived above. If the small 
grains accumulate, only of order 1 -MEarth of crystalline 
material has to be set free in collisions, a number we can 
comfortably account for. 
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We note that the crystalline dust in zone 2 of the small 
grain component is distributed in the model fit over a 
larger range than the amorphous dust. One way in which 
this could be realized is that small dust grains produced 
by the collisional cascade will be driven out by radia- 
tion pressure. Because the disk is flared, the grains will 
not be able to leave the system, but will "land" further 
out on the disk surface. There they will enrich the super- 
heated disk surface layer which is known to be responsible 
for the features in the spec tra of irradiated flaring disks 
llChiang fc Goldreichl [l997h . This redistributi on of pro- 
cessed dust is similar to the model proposed bv lShu et alJ 
( 199<|) for the origin of Chondrules and calcium- aluminum 
rich inclusions (CAIs) in meteorites. In the IShu et ail 
model, grains are lifted off the disk by the jet which is leav- 
ing the disk very close to the star in the active accretion 
phase of the disk. The grains then are heated up to melt- 
ing temperatures by the stellar light and transported to 
larger distances by radiation pressure. In what we propose, 
the grains are produced by fragmentation of larger bodies, 
and are scattered out of the disk. Small grains for which 
radiation pressure is dynamically important will be driven 
away from the star and by a combination of gravity and 
disk geometry land further out on the disk surface. This is 
the location where we see these grains today. A schematic 
representation of the above sketched disk structure and 
processes is plotted in Fig. ^3 

An alternative scenario for the in situ formation of 
the crystalline s i licate s is annealing by nebular shocks. 
lHarker fc Descbl l)2002|) argue that local gravitational in- 
stabilities, likely to exist in protoplanetary disk beyond 5 
to 10 AU, could cause shocks. These shocks could suffi- 
ciently heat small (< lpm) amorphous grains in order to 
crystallize them. Interestingly, a developing planet, open- 
ing a gap in the disk, would produce sh ock waves through 
tidal intera ction with the disk (e.g. iLin fc Papaloizoul 
Il980lll993h . However, the thus generated shocks should 
radially extend over the very wide observed region of the 
forsterite grains. It remains to be seen whether the gas 
densities are high enough in regions extending far beyond 
10 AU for this mechanism to be effective. 

The consequence of the above proposed scenario is that 
Hale-Bopp, given the very similar dust composition, may 
contain material formed by collisions of bodies formed ear- 
lier on in the history of the Solar system. The thermal con- 
tact of the dust species in Hale-Bopp implies that we are 
not witnessing the destruction of Hale-Bopp type objects 
in the HD 100546 system, but the creation of the material 
similar to that which formed Hale-Bopp before the dust 
grains coagulated. Interestingly this would place the for- 
mation of Hale-Bopp after the formation of one or more 
of the gas planets in our solar system. 

6. Summary 

We can summarize the results of our analysis of the ISO- 
spectra of HD 100546 and comet Hale-Bopp as follows: 



1. The ISO/SWS spectra of HD 100546 and comet Hale- 
Bopp can be fitted with very similar dust composi- 
tions and properties. In both cases, the bulk of the 
material is dominated by amorphous silicates with car- 
bonaceous grains as the second most important con- 
stituent. In addition, significant amounts of forsterite, 
crystalline water ice, metallic iron and iron sulfide are 
present. 

2. There is one major difference in the dust proper- 
ties of the two investigated objects. The chemical 
constituents of the dust in Hale-Bopp are found to 
be in thermal contact while this is not the case for 
HD 100546. 

3. The fraction of crystalline silicates, i.e. forsterite 
and/or very iron poor olivine, in HD 100546 increases 
with decreasing temperature, i.e. with increasing ra- 
dial distance from the central star. Such a distribution 
is not in agreement with predictions by radial mix- 
ing models. In such models the crystals are formed by 
thermal annealing in the very inner parts of the disk 
and are subsequently transported out and mixed with 
amorphous material. 

4. HD 100546 shows a component of small grains (< 10 
/xm) with temperatures of about 200 K that emits in 
the mid-IR and that is not present in the spectra of 
other Herbig Ae/Be stars. 

5. At about the distance where the bulk dust temperature 
is ~ 200 K, i.e. at ~ 10 AU, the vertical extent of the 
disk must be in excess of that expected in a standard 
flaring disk model. This is required to explain the large 
fraction of stellar light intercepted and re-emitted by 
grains emitting in the mid-IR. 

6. On the basis of the above three results, we propose the 
following scenario: the crystalline dust is produced by 
collisional destruction of differentiated objects. These 
collisions are probably induced by gravitational in- 
teraction with a proto- Jupiter. The grain component 
emitting at ~ 200 K may be identified by the dust ex- 
pected to be produced in the collisional cascade that is 
thought to occur when large bodies collide. This pre- 
dominantly crystalline material may be driven out by 
radiation pressure and "land" further out on the (flar- 
ing) disk surface. The gap cleared out by the proto- 
Jupitcr allows direct stellar light to produce an ex- 
tended rim or wall at the far side of the gap. It is 
thought to be this rim that intercepts the radiation 
responsible for the strong mid-IR component in the 
spectrum of HD 100546. 

7. The similarity in spectral appearance between 
HD 100546 and Hale-Bopp implies that similar pro- 
cesses in the early history of our solar system may 
have been relevant for the formation of Hale-Bopp. The 
comet may have assembled from debris dust produced 
through collisions of differentiated objects. 
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